Splashing of liquids: interplay of surface roughness with surrounding gas 
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We investigate the interplay between substrate roughness and surrounding gas pressure in con- 
trolling the dynamics of splashing when a liquid drop hits a dry solid surface. We associate two 
distinct forms of splashing with each of these control parameters: prompt splashing is due to surface 
roughness and corona splashing is due to instabilities produced by the surrounding gas. The size 
distribution of ejected droplets reveals the length scales of the underlying droplet-creation process 
in both cases. 

PACS numbers: 47.55.Dz, 47.20.Cq, 47.20.Ma 



I. Two contributions to splashing 

A liquid drop hitting a surface will often splash, rup- 
turing into many tiny droplets. Although splashingis 
ubiquitous and important in many applications [J, [Hibl , 
the control parameters governing its occurrence have not 
yet been fully explored. Impact velocity, surface tension, 
viscosity and substrate roughness have long been known 
to be important 0, S @| • However only recently was it 
discovered that splashing can be completely eliminated 
on smooth surfaces simply by lowering the surrounding 
gas pressure In this paper we focus on the role that 
surface roughness plays in producing a splash. Given 
the discovery that the gas pressure can, by itself, cause 
splashing, it is apparent that only by removing the sur- 
rounding gas can the other control parameters, such as 
surface roughness, be investigated in isolation. In order 
to make progress, we therefore eliminate all vestiges 
of the corona splash by working at low gas pressures. 
This insures that the effects of surface roughness are 
unperturbed by the effects of the surrounding gas. 

Two types of splashing are well-known in the litera- 
ture Q: "corona" splashing, where a corona forms and 
subsequently ruptures and "prompt" splashing, where 
droplets emerge at the advancing liquid-substrate con- 
tact line. Although the two phenomena are distinct there 
has been no fundamental understanding of their separate 
causes. The results reported here suggest a simple expla- 
nation: corona splashing is due to the presence of sur- 
rounding gas and occurs above a critical gas pressure and 
prompt splashing occurs on rough surfaces even in the ab- 
sence of air. By studying the ejected-droplet size distri- 
bution of these two forms of splashing, we can determine 
the characteristic length scale of the droplet-creation pro- 
cess in each case. This corroborates our interpretation of 
how the splash is formed. 

Our previous experiment 0] clearly showed that sur- 
rounding gas causes the corona splash on smooth sur- 
faces. Here we study the effect of another control pa- 
rameter: surface roughness. In order to do this cleanly, 
we must lower the gas pressure until the effects of air 
are negligible. We achieved this by doing an experiment 
in a helium atmosphere at 13kPa pressure, much lower 



than the threshold pressure [7[ for the surrounding gas 
to produce a corona splash. In this situation, therefore, 
splashing is caused entirely by surface roughness. 

In our experiments, we use ethanol which has den- 
sity p = 0.789g/cm 3 , viscosity v = 1.36c5i and sur- 
face tension a — 22mN/m. We filmed drops, of di- 
ameter D = 3.4 ± 0.1 mm, released inside a transpar- 
ent vacuum chamber from a nozzle at a height 95 cm 
above a rough substrate. The impact speed of the drop, 
Vq — 4.3±0.1m/s, was determined by analyzing the drop 
position in subsequent frames of each movie. Rough sub- 
strates were obtained by using high-quality sandpapers 
uniformly coated with closely packed particles (microcut- 
paper-discs® , Buehler Ltd.). The roughness, R a , defined 
as the average diameter of the particles, was varied be- 
tween 3/zm and 78/im. Clean glass microscope slides were 
used as the smooth substrates. The gas pressure, P, was 
varied. Ethanol wets our substrates and does not re- 
bound after hitting the surface [?|. We also note that 
although ethanol wets the substrate (zero degree static 
contact angle), the profile of the liquid film during ex- 
pansion has a finite thickness at the edge, instead of de- 
creasing smoothly to zero. All experiments were done at 
constant temperature, 23.6 °C 

As we increase the roughness, i? Q , we see an evolution 
in behavior as shown in Fig. 1. We compare the splash 
created solely by surface roughness at low pressure with 
the splash created when air is also present. In each panel 
the top row is at low pressure, P = ISkPa, and the bot- 
tom row is at atmospheric pressure, P — lOOkPa. Fig. la 
shows the result for a smooth surface: at low pressure 
there is only a liquid film expanding smoothly on the sub- 
strate with no splashing whereas at atmospheric pressure 
there is a corona that breaks up into many small droplets. 
This shows unambiguously that air causes corona splash- 
ing. Fig. lb shows the result for a small amount of 
roughness, R a = h[im. In this case, at low P there are 
two regimes: an early stage with prompt splashing which 
is followed by a peaceful regime where no splashing oc- 
curs. At atmospheric pressure, there is a single regime 
throughout the expansion that resembles corona splash- 
ing except that it is not as symmetric as when the surface 
is completely smooth. Thus, for roughness R a = 5fj,m , 
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FIG. 1: Photographs of splashing as a function of gas pressure 
and surface roughness. The left and right columns are 0.2ms 
and 0.5ms from the time of impact. For each value of surface 
roughness, the top panel is at a low pressure, P — 13kPa and 
the bottom panel is at atmospheric pressure, P = 100/cPa. 
(a) Splash on a smooth substrate which is a clean microscope 
slide, (b) Splash on a substrate with roughness R a — 5fim. 
(c) Splash on a substrate with roughness R a = 78p,m. 



there is a clear difference as P is increased, showing that 
surrounding gas can still be important at atmospheric 
pressure. For large roughness, R a — 78/ito, we see typ- 
ical prompt splashing, where droplets are ejected from 
the expanding contact line during the entire expanding 
process for both low and atmospheric pressures. 

In summary, there are two contributions to splashing: 
a corona part caused by gas and a prompt part caused by 
surface roughness. At low pressures, for small roughness 
(R a = 5/xm), there is no corona and there is only a small 
amount of prompt splashing produced at the beginning. 
This prompt splashing disappears at later times. For 
large roughness (R a — 78pm) , splashing at the advancing 
contact line is produced throughout the entire duration of 
the expansion. At atmospheric pressure there is a tran- 



sition as the roughness is increased: the corona splash 
dominates at small surface roughness and the prompt 
splash dominates at large R a . 

These results indicate that in the absence of air, splash- 
ing is caused when the expanding liquid film, of thickness 
d, becomes destabilized by surface roughness; but if the 
roughness is too small or the film is too thick then no 
splashing will occur. Initially the expanding film thick- 
ness d is of molecular size and increases in thickness 
during expansion as liquid is added to the film. Thus 
for small roughness (R a = bum), splashing occurs im- 
mediately after impact and is followed by a quiescent 
stage as the film becomes much thicker than the rough- 
ness. For large roughness (78pm), splashing continues 
throughout the film expansion since d never grows large 
enough to be unperturbed by the roughness. From the 
photographs, we estimated the liquid film thickness, d, 
at the point where splashing stops for R a = bpm and 
found d ~ 50/iTO. This suggests the following criterion 
for prompt splashing: 

^ = C(We,Re) (1) 

where C(We, Re) is a dimcnsionlcss number depending 
on Weber number We = pV^D/a and Reynolds number 
Re = pVoD/p. F° r the impact conditions in Fig.l, 
We « 2400 and Re w 11500, we conclude C ~ 0.1. Fur- 
ther studies are necessary to establish the dependence 
of C on We and Re. For example, at sufficiently low 
velocity we might not expect to see any splash at all. 
We emphasize that C(We,Re) can only be measured 
accurately when the effect of gas pressure is negligible 
since otherwise the corona component will contaminate 
the results. This shows the importance of separating the 
two types of splashing by working at small P. 

II. Size distribution of prompt splashing 

After impact, the liquid breaks up into many tiny 
droplets. The size distribution of these emitted droplets, 
N(r), may retain an imprint of the droplet creation 
process and thus provide a clue to the mechanism 
initiating the interfacial instability. Here we report N(r) 
for prompt splashing on rough substrates. 

We measure the size of emitted droplets by adding a 
small amount of ink (Sanford® black stamp pad inker) 
to the ethanol, at the volume ratio 1:6 (ink:ethanol). The 
mixed liquid has the following material properties: den- 
sity: p = 0.833±0.002#/cm 3 , viscosity: v = 3A±0.2cSt, 
surface tension: a = 23 ± 3mN/m. (a is measured both 
in air and helium atmosphere under low pressure. Both 
cases give the same result, 23 ± 2>mN/m.) Except for 
the viscosity, which increases by a factor of 2.5, these 
values are close to those for pure ethanol. We have also 
checked repeatedly with high-speed video and found that 
the splashing pictures look very similar with and with- 
out ink, therefore we made sure that the addition of ink 
does not change splashing qualitatively. We surround the 
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FIG. 2: Calibration curve of stain area as function of droplet 
volume. Error bars come from the fluctuation of stain areas 
for the same size droplet. The fitting curve is: y = 233 * 
x - 67 — 119. y is the stain area in units of pixels, and x is the 
volume of the drop in units of nanoliters. 



impact point with a cylinder of diameter 8.89 ± 0.02cm 
rolled from a sheet of white paper. After splashing, the 
droplets hit the cylinder, leaving stains. 

With careful calibration, we convert the sizes of the ink 
spots on the paper to the sizes of the ejected droplets. 
Our calibrations work well for drops with radius larger 
than r = 50/im. Below that radius, our resolution is in- 
adequate to obtain a reliable calibration. We used a 0.5/xl 
syringe to deliver tiny droplets of the ethanol/ink fluid 
with known volume to white paper. By measuring the 
area of the stains created by these droplets we obtained 
the calibration curve in Fig. 2. We checked the effect of 
droplet velocity on the stain area and found no significant 
effect within our experimental accuracy. By checking the 
shape of a stain, we can ascertain whether that spot was 
caused by a single ejectile or by two separate drops that 
landed in overlapping locations. We found less than 4% 
overlapping stains. We excluded them from our distribu- 
tion curves. 

We first determine N(r) for prompt splashing on rough 
substrates. To ensure that these distributions are due 
solely to prompt splashing, we performed the measure- 
ments at low pressure. Fig. 3a shows N(r) for several 
values of roughness, R a . To obtain good statistics, each 
distribution is an average over ten to twenty experiments 
taken under the same conditions. We added the number 
of spots found in all the experiments and then divided the 
total counts by the number of experiments. The straight 
lines in the figure indicate that N(r) decays exponentially 
with a characteristic decay length, r : 



N ~ exp(—r/ro) 



(2) 
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Fig. 3b shows that the decay length of these lines, ro, 
increases with increasing R a . At small R a , the decay 



FIG. 3: The distribution of ejected drops in a prompt splash 
on rough substrates, (a) N(r) versus r for splashes on sub- 
strates with three values of roughness: R a — 16/im(«), 
R a = 5fim(o), and R a = 3/im(A). The exponential fitting 
functions are: — , exp(-r/0.023); — , exp(-r/0.006); and 
— , exp(-r/0.004). (b) The decay constant, ro, of the expo- 
nential decay in N(r), as a function of substrate roughness 
R a . For small values of roughness, the decay constant is ap- 
proximately linear in the roughness. At large roughness, the 
decay constant saturates. The sizes of particles on sandpa- 
per are randomly distributed around an average value, R a - 
The fluctuation of particle sizes gives us R a error bars; and 
the standard deviation from exponential fit gives the ro error 
bars. Inset of (b) shows one sheet of paper with ink spots pro- 
duced by R a = 25/im. Most spots are randomly distributed 
within a narrow band. 



length is close to the value of the roughness: ro ~ R a - 
At large roughness, this relationship breaks down as ro 
appears to saturate at a constant value. 

This behavior is consistent with the prompt splash- 
ing criteria, Eq. 1. When the roughness is small, the 
drop stops splashing at a film thickness determined by 
R a - This sets the correlation between the droplet de- 
cay length and roughness: ro ~ R a - However, when the 
roughness is too large, the drop never stops splashing, 
and ro can only increase up to the maximum thickness 
of the expanding liquid film. From the photographs we 
estimate thickness d ~ 100/im. Fig. 3b shows that the 
saturation occurs at approximately 40/xm, which is about 
the same order of magnitude as the thickness of the liquid 
film at its terminal position. 

How robust is the decay length, ro, with respect 
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FIG. 4: Size distribution of different impact velocities. For 
a fixed roughness R a = 25/im, two impact velocities were 
tested: V = 4.3 ± 0.1m/s(o) and V = 5.2 ± 0.1m/s(i). (a) 
shows N(r) in linear4inear scale. It's clear that the higher 
velocity case produces about 20% more splash, (b) plots the 
same data in log-linear format. The two data sets now seem 
much closer only because of the log4inear way we plot them. 
We can fit both curves with the same functional form: A ■ 
exp(— r/0.025)(the solid line), by only varying the amplitude, 
A. This implies that ro = 0.025mm is independent of impact 
velocity. 



to variations of the impact velocity Vq, drop size 
R, or surface tension a? In Fig. 4 we show the 
size distribution for two different impact velocities: 
V Q = 4.3 ± 0.lm/s(We = 2155, Re = 3648) and 
V = 5.2 ± 0.1m/ s(We = 3151, Re = 4412). Fig. 
4(a) shows that there is significantly more splashing 
at higher impact velocity. Fig. 4(b) plots the same 
data in log-linear format and shows that both curves 
can be fitted by straight lines with the same decay 
length, r . This implies that r is independent of Vq. 
Further experiments addressing other parameters should 
be done in the future. Although the absolute value 
of r might change with those parameters, we expect 
that the two general features: (1) the distribution 
decays exponentially and (2) ro increases with surface 
roughness, would still be valid. 

III. Size distribution of corona splashing 




FIG. 5: Size distribution of ejected droplets in a corona splash 
at high pressure. Upper inset shows the corona film before 
it breaks up. From this picture we can estimate the film 
thickness to be 20 ~ 40pim by measuring the thickness of 
the corona rim. (Although rim should be thicker than the 
film, we assume that the film and the rim are approximately 
the same size.) Lower inset is a reproduction of the sheet of 
paper with the ink spots showing that the ejected droplets hit 
the paper at random locations over a large area. Main panel 
shows the number of droplets of a given size per impact, N(r), 
versus droplet radius, r, for a corona splash at two pressures: 
P = lOOkPa (•) and P = 80kPa (x). The exponential fitting 
functions are respectively: — , exp(-r/0.020); and — , exp(- 
r/0.014). 



We have also performed measurements of the size dis- 
tributions produced in a corona splash. This allows us to 
compare prompt and corona splashing. Fig. 5 shows the 
results of splashing on a smooth surface at two pressures, 
lOOkPa and 80kPa, both in the high pressure regime. The 
inset shows the spots created by the droplets ejected from 
a single splash. The spots are randomly distributed. The 
main panel shows the distribution of droplet sizes, N(r), 
for the two pressures. We again see an exponential dis- 
tribution at large radius r, N ~ exp(— r/ro), indicating 
the existence of a characteristic length scale. The values 
we find for r (given in the caption of Fig. 5) arc com- 
parable to the corona thickness which we estimated from 
the movies to be between 20 and 40/im. (The corona 
thickness was estimated by measuring the edge of the 
corona in pictures similar to the one shown in the inset 
to Fig. 5.) This suggests that the corona thickness de- 
termines r . That ro decreases as the pressure decreases 
suggests that lower pressure leads to less splashing and a 
thinner corona. The data in Fig. 5 shows that the size- 
distribution data can probe the slight differences in the 
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FIG. 6: A splash just above the threshold pressure, Pt- (a) 
Images of a drop splashing show the ejection of droplets and 
the undulations in the expanding rim. The times shown are 
measured with respect to the time of impact, (b) Inset is a 
reproduction of the sheet of paper showing that the droplets 
hit the paper in a well-defined horizontal line. Main panel 
of (b) shows the Fourier transform of the lateral positions 
of the spots in the inset. The peak indicates a well-defined 
spacing between the ejected droplets at Pt- (c) The number of 
droplets at fixed sizes, N(r), is plotted versus droplet radius, 
r, at pressures slightly higher than Pt- P 34kPa(«), P= 37 
kPa(o) and P = 39kPa(x). The solid line is a Gaussian fit 
centered at ro = 0.11mm. The peak in N(r) shows that the 
average size of the droplets is approximately the size of the 
undulations at the rim of the expanding film. 



corona formed at the two pressures. 

Both prompt and corona splashing have exponential 
distributions for N(r). Although there is a characteristic 
length scale in both cases, the control parameters 



governing the length scales are different: roughness in 
Fig. 3 and gas pressure in Fig. 5. However, at the 
threshold pressure [3] for the corona splash on a smooth 
surface, there is a qualitatively different distribution 
of droplet sizes. Fig. 6a shows there is no corona and 
that discrete droplets emerge from the expanding liquid 
which has periodic undulations along its rim. The ink 
spots shown in the inset at the lower left corner of Fig. 
6b have the striking feature that they fall in a horizontal 
line, indicating that the ejected droplets have the same 
angle between their trajectories and the substrate. In 
addition, the spots are approximately equally spaced. 
The main panel of Fig. 6b shows that the Fourier 
transform of the lateral positions of the spots has a peak 
indicating this spatial order. Moreover the spot sizes are 
more uniform than those seen at high pressures in the 
inset to Fig. 5. Fig. 6c shows N(r) at pressures close 
to Pt- The peak at ro = 0.11 mm indicates that most 
droplets are about the size of the rim undulations. 

IV. Conclusions 

By controlling the gas pressure and surface roughness 
we have identified two mechanisms for splashing: the sur- 
rounding gas is responsible for corona splashing and sur- 
face roughness is responsible for prompt splashing. This 
explains the long-standing puzzle about why two distinct 
types of splashes exist. We have also found that there are 
characteristic lengths in the distribution of ejected drops 
in both cases. Similar exponential dependence is found 
in ligament breakup, implying a possible connection with 
that process[l(| EJ. However, the exponential depen- 
dence of droplet sizes found in the splashing of liquids is 
in contrast to what is found in the shattering of a solid. 
In that case, there is a power-law distribution of sizes 
of the shattered fragments and no characteristic length 
scale [HI, [l3j]. The characteristic lengths we find in liq- 
uid splashing experiments reveal the microscopic length 
scales associated with the droplet creation process. On 
rough surfaces, the length scales we determine are con- 
sistent with our interpretation of the prompt splash in 
terms of the surface roughness. These results also sug- 
gest a means for controlling the sizes of ejected droplets 
in a splash. 
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